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ABSTRACT 

Long  waves  are  affected  by  bottom  topography  and  under  certain 
conditions  may  be  trapped  along  topographical  contours  which  then  act  as 
wave  guides  transmitting  wave  energy  for  great  distances  with  little  loss. 
This  study  examines  waves  trapped  along  a  submerged  ridge  described  by 
straight  parallel  bottom  contours  which  in  cross  section  are  composed  of 
constant  slope  segments  bounded  on  either  side  by  constant  depth  segments. 
Solutions  are  found  for  time  harmonic  waves  periodic  in  the  along  ridge 
direction  and  of  exponential  decay  behavior  normal  to  the  ridge  over  the 
constant  depth  segments.  Over  the  linearly  varying  topography  describing 
the  ridge,  the  solution  is  in  terms  of  two  Kummer  (or  Whittaker) 
functions.  For  a  given  geometry,  a  dispersion  equation  is  obtained 
relating  the  wave  frequency  to  the  along  ridge  wave  number  for  trapped 
waves.  A  constant  Coriolis  parameter  is  included,  but  primary  interest  is 
on  class  I  (high  frequency)  waves.  A  comparison  of  cutoff  frequencies 
predicted  for  this  piecewise  continuous  ridge  and  those  for  a  segmented 
constant  depth  ridge  is  made,  and  the  appropriate  scaling  factors  between 
the  two  results  are  discussed. 
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INTRODUCTION 

The  study  of  topographically  trapped  waves  appears  to  have  begun 
with  Stokes’  description  of  edge  waves  on  a  linear  topography  of  straight 
parallel  bottom  contours,  Stokes  (1846).  While  this  solution  was 
originally  felt  to  be  of  only  theoretical  interest,  observations  during 
the  last  few  decades  have  given  new  physical  significance  to  these  forms 
of  waves.  This  study  shall  examine  the  influence  of  one  dimensional 
(parallel  contours)  bottom  topography  on  the  trapping  and  guiding  of  ocean 
wave  energy.  In  particular,  linearized  long  wave  theory,  for  a 
homogeneous  perfect  fluid,  with  rotation  is  used  over  topographies  which 
consist  of  constant  and  linear  depth  cross  sectional  variations. 

Before  further  detail  on  this  study,  a  review  of  some  background  material 
is  in  order.  LeBlond  and  Mysak  (1978)  provide  an  important  overview  of 
the  area  of  ocean  waves;  more  specific  results  are  given  below. 

Eckert  (1951)  showed  that,  for  shallow  water  theory,  the  Stokes  solution 
was  just  the  lowest  mode  of  an  infinity  of  modes  and  Ursell  (1952)  obtained 
these  results  generalized  to  three  dimensional  linearized  water  wave  theory 
leading  to  a  mixed  spectrum  of  discrete  and  continuous  eigenf requencies . 
Reid  (1958)  showed  that  the  earth’s  rotation  would  have  a  modifying  effect 
on  these  class  I  or  inertio-gravitational  waves  (which  can  propagate  in 
both  directions  along  a  N-S  coast  albeit  with  differing  phase  speeds  for 
non-zero  rotation)  and  would  also  allow  for  a  class  II  or  quasi-geostrophic 
wave  which  can  propagate  in  only  one  direction  and  can  not  exist  in  the 
absence  of  rotation. 

The  differences  between  class  I  (high  frequency)  and  class  II  (low 
frequency)  shallow  water  waves  has  lead  to  two  roughly  distinct  classes  of 
literature — those  primarily  concerned  with  relatively  nearshore  phenomena, 
e.g.,  tsunami  periods  and  shorter  (under  1  hour),  see  Munk,  Snodgrass  and 
Carrier  (1956),  and  those  concerned  with  oceanographic  scale  phenomena, 
e.g.,  periods  of  the  order  of  days.  Clearly,  these  interests  overlap  but 
a  distinction  can  be  made  on  the  basis  of  the  wave  frequency,  w,  versus 
the  Coriolis  frequency,  f.  If  w  >>  f,  it  is  the  "edge”  wave  or  class  I 
wave  which  is  usually  of  primary  interest  while  if  u  «  f,  it  is  the 
"quasi-geostrophic"  or  class  II  wave.  This  study  is  directed  at  the  former; 
the  second  class  of  waves  will  be  discussed  in  a  later  paper. 


It  is  clear  from  the  various  formulations  of  these  problems  that 
waves  may  be  trapped  over  completely  submerged  topography,  as  shown  by 
Hidaka  (1976) — see  Meyer  (1971)  for  a  review  of  resonance  of  unbounded 
water  bodies,  Shen,  Meyer  and  Keller  (1968)  for  a  geometrical  optics 
approach  and  Buchwald  (1968)  for  a  solution  of  the  piecewise  constant  depth 
case  and  a  discussion  of  the  general  nature  of  such  problems  as  solutions 
to  Sturra-Liouville  equations.  One  aim  of  the  present  study  is  to  see  how 
well  the  piecewise  constant  depth  solution  represents  solutions  for  gently 
sloping  bottom  topographies  by  comparison  to  another  analytical  solution 
for  piecewise  linear  depths*  The  solution  to  this  linear  topography  was 
noted  by  Shaw  (1974)  and  Guza  and  Davis  (1974)  in  terms  of  Rummer  functions, 
but  no  numerical  results  were  given  until  Shaw  (1977),  due  to  an  apparent 
lack  of  tables  and/or  subroutines  for  these  functions. 

While  the  question  of  resonances  is  difficult  enough  for  the 
trapped  wave  problem,  the  prior  question  of  how  energy  arrives  to  be 
trapped  must  also  be  considered.  For  trapped  waves,  the  energy  may  be  supplied 
by  nonlinear  effects,  local  atmospheric  forcing  or  even  local  seismic 
disturbances  since  a  linearized,  unforced  theory  only  provides  for  the 
persistence  of  such  trapped  motions  once  begun.  However,  it  must  be  pointed 
out  that  total  wave  trapping  is  exceptional;  Longuet-Higgins  (1969)  has  shown 
that  curved  topographies  lead  to  ,f leaky”  trapped  modes.  Rather  than 
negating  the  importance  of  trapped  wave  studies,  these  results  enhance  them, 
since  this  provides  a  simple  mechanism  for  these  cases  by  which  incident 
wave  energy  can  be  fed  into  the  resonating  system.  Alternatively,  for  a 
ridge  of  finite  length,  wave  energy  may  enter  the  "system"  at  the  ends, 
exciting  a  resonance  along  the  ridge  as  well  as  across  the  ridge. 

This  mechanism  may  be  particularly  appropriate  at  locations  where  a  long 
ridge  terminates  at  a  coastline,  e.g.  the  Chattham  rise  off  of  the  east  coast 
of  New  Zealand,  Heath  (1979). 


FORMULATION 


The  problems  considered  here  are  concerned  with  wave  energy 
trapping  by  a  topography  with  either  a  one  dimensional  piecewise  discon¬ 
tinuous  cross  section  composed  of  constant  depth  segments  or  a  piecewise 
continuous  cross  section  composed  of  linear  depth  segments •  Although 
rotation  of  an  f  plane  will  be  included,  emphasis  will  be  on  first  class 
surface  gravity  waves  on  a  homogeneous  fluid  where  Coriolis  effects  are  a 
modifying  rather  than  a  fundamental  influence. 

The  vertically  integrated  equations  of  motion  in  a  coordinate 
system  with  contours  parallel  to  the  y  axis  are 


a) 

3u/3t  -  fv  =  -g  3£/8x 

b) 

8v/3t  +  fu  =  -g  8£/8y 

(i) 

c) 

8 [Hu] / 8x  +  8[Hv]/8y  =  -8c/8t 

with  water  depth  H  and  Coriolis  parameter  f. 

Assuming  time  harmonic  behavior,  exp(-iu)t),  and  periodic  along- 
contour  behavior,  exp(iky),  allows  these  equations  to  be  solved  separately 
for  the  horizontal  velocities  u,  v  and  free  surface  elevation, 

a)  u  =*  [-iwg  8c/8x  +  ikgc]/[u>2  -  f2] 

b)  v  =  [uigkc  +  fg  3£/3x]/[u)2  -  f2]  (2) 


c) 


a  SB  -  Hk* 

a)  dx 


C  =  0 


These  equations  may  be  non-dimensionalized  with  respect  to  some  reference 
depth  and  length  L,  i.e.,  using 

?  =  C/%,  h  =  h/hr 

(x,y)  =  (x,y)/L,  k  =  kL 
(u,v)  -  (u,v)  ^wL/gHj^j 
f  -  f/w,  2  »  wL/^g^j1/2 


equations  (2)  become 


a  -  i[fkc  -  3c/3x]/[l  -  f2] 
V  -  [ki  +  f 3c/35c]/ [1  -  f2] 

Km  +  dK<*>  f!i  +  r 


k?  ^ 
dx 


hk2 


C  “  0 


which  will  have  different  solutions  for  different  topographies,  h(x). 

The  simplest  case  arises  for  constant  depth,  e.g.,  Hq  (not  necessarily  the 
reference  depth)  such  that  hQ  =  Hq/Hr, 

5o  +  [*2<1  -  I2)  -  Koi‘2)«  '  0 


leading  to 

£  *  A  exp  [+Xx]  +  B  exp[-Xx] 

where  X  is  defined  by  the  indicial  equation 
X2  -  k2  -  u)2(l  -  f2)/h 

o 

For  values  of  and  f  such  that  k2  >  532  (1  -  f2)/h^,  X  is  real  and  the 
solution  for  £  is  exponential  in  x  while  for  k2  <  32(1  -  f2)/hQ,  X  is 
imaginary  and  the  solution  is  sinusoidal  in  x. 

If  a  linear  depth  is  considered,  h(J5)  *  [6  +  yx]  ,  this  equation 


reduces  to 


I5  +  yx]  ^=7  +  Y  ^4  +  [0)2(1  -  f2)  -  kfy  -  (<S  +  yx)k2]c 


Introduce  z  -  6  +  yx;  this  equation  is  then 
d2£  ,  d£  ,  fo)2(l  -  f 2 )  -  kfy  k2 


2f|  + 

dz * 


d£  .  fs2(l 

dz 


“  T  2 


'  l  -  0 


which  is  reducible  to  the  confluent  hypergeometric  equation,  e.g.,  Erdelyi, 
et  al.  (1953).  Solutions  to  this  equation  are  of  the  form,  with  £  >  0: 


£  -  expf-£/2] 


{ClM(a,l,0  +  c2U(a,l,ol 


where  £  *  -2kz/y,  a  =  (1  -  f ) / 2  +  u)2(l  -  f2)/2yk,  (y  ^  0)  and  M  and  U  are 
Kummer  Functions,  with  arbitrary  coefficients  c^  and 

For  £  <  0,  the  Kummer  transformation  leads  to  an  alternate  form: 

l  -  exp[+£/2]jc’  M(a\l,-0  +  c'2  0(aM,-()j 

where  a’  *l-a~i  +  ~-  -■  Z—  — 

2  2  2yk 

A  program  for  Kummer  functions  of  both  kinds  is  given  by  Shaw  and 
Neu  (1977)  in  an  internal  SUUY-B  report  and  is  available  from  the  authors. 
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STEP  RIDGE  TOPOGRAPHY  -  DISCONTINUOUS,  CONSTANT  DEPTH  CROSS  SECTION 

This  problem  was  discussed  by  Buchwald  (1968)  and  forms  a  test 
case  against  which  the  later  continuous  depth  case  may  be  examined. 

Since  constant  depth,  discontinuous  (piecewise  continuous)  topographies 
are  much  simpler  to  use  than  any  other  topography,  it  is  of  great  practical 
interest  to  know,  analytically,  how  well  they  may  represent  a  gradually 
varying  topography.  The  basic  comparison  will  be  done  for  the  dispersion 
equation,  in  particular  for  the  low  frequency  cutoff  frequencies  and  other 
significant  properties. 

Consider  a  cross  section  defined  by  three  regions  of  constant  depth 

Region  I:  -«  <  x  <  -B;  H(x)  *  H^ 

Region  II:  -B  <  x  <  +B;  H(x)  *  H^ 

Region  III:  +B  <  x  <  +»;  H(x)  *  H^ 

where  coordinates  are  chosen  such  that  H^  £  H^. 

Using  B  and  H^  as  reference  lengths  for  horizontal  and  vertical 
distances  respectively,  the  solution  in  Region  I  is 


-  Ai  exp[xi  s] 


where  X  *  [k2  -  w2(l  -  f2)]1/2  and  the  portion  of  the  solution  which 
diverges  as  x  is  surpressed.  In  Region  II,  an  oscillatory  behavior  is 


anticipated,  i.e. 


*  A2  cos 


+  B^  sin 


M 


where  *2  *  |w2(l  -  f2)/h^  -  k2J*/2.  Finally,  in  Region  III  we  anticipate 
another  exponentially  decaying  solution 


4ll  "  A3  exp[  X3X] 
where  X^  «  ^k2  -  w2(l  -  f2)/h^1/2 . 

Clearly  k  is  bounded  by  aj£(l  -  f2)/!^1/2  from  above  and  by 
w(l  -  f2]1/2  from  below. 


Continuity  of  surface  elevation  and  mass  flow  requires 


Cj  (x  =  -1) 

II 

hi  =  _1^ 

Uj.  (x  =  -1) 

=  H2 

°n  (x  =  ’1} 

Cj-jCx  =  +1) 

88 

Cm(x  =  +1) 

h2  u^Cx  =  +i) 

=  R3 

5III(*  =  +1) 

leading  to  a  system  of  4  homogeneous  linear  algebraic  equations  on 

Ja^,  9  ,  A^j .  The  determinant  of  the  coefficient  matrix  of  this  system 

must  vanish  in  order  to  have  a  non-trivial  solution;  this  defines  the 
dispersion  equation.  For  simplicity,  consider  the  symmetric  case,  =  H^, 
without  rotation,  f  =  0.  The  dispersion  equation  then  reduces  to 


2h 2  cos  ^2  =  (^2  *2  ~  ^l)Sin  2^2 

Along  the  line  c3  =  k  which  forms  a  lower  bound  on  the  allowable 
values  of  k  for  trapped  waves,  =  0  and  the  dispersion  equation  reduces 
to 

sin  2X^  =  0,  *2  ~  n7T^*  n  =  0»1»2,3.  .  . 

These  are  the  cutoff  points,  corresponding  to 


=  (nir /2)  jh 2/  ^l-h2j  |  ^ 
*  (nn / 2 ) |h2/ ^1— h 2]  j  1 


The  shape  of  the  free  surface  separates  into  symmetric  and 
antisymmetric  modes.  For  even  n,  is  zero  and  =  A3  8  cos (nn/ 2) A^ , 

i.e.  these  are  the  symmetric  modes,  while  for  odd  n,  A^  is  zero  and 
A^  *  -A^  ■  -sin(mi/2)A2>  i.e.  the  antisymmetric  modes. 
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LINEAR  RIDGE  TOPOGRAPHY  -  CONTINUOUS  LINEAR  CROSS  SECTION 

Consider  a  ridge  composed  of  a  piecewise  continuous,  linear  cross 
section  defined  by 

Region  I:  <  x  <  -A;  H(x)  * 

Region  II:  -A  <  x  <  0;  H(x)  =  H2  -  ^  -  H2jx/A 

Region  III:  0  <  x  <  B;  H(x)  =  H2  +  ^  -  H2jx/B 

Region  IV:  B  <  x  <  ®;  H(x)  =  H3 

as  shown  in  Figure  1.  A  suitable  reference  depth  is  and 

reference  length  L  -  A. 

Then  in  Region  I,  h^  =  1  and 

=  A.^  expj^Jtj  +  B.^  exp  £-Ajxj 

where 

\1  =  [k2  -  S2(l  -  f2)]1/2 

Clearly,  B^  must  be  set  equal  to  zero  and  k2  >  ui2(l  -  f2)  to  have  an 
exponentially  decaying  wave  as  x  •+  -*>. 

Similarly  in  Region  IV,  h  ■  H^/H^  =*  h^  and 

4v  =  A4  exP  [_X4x]  +  B4  exP  [+X4x] 

where 

A4  =  [k2  -  m2(l  -  f2)/h3j1/2 

Again  we  require  a  decaying  solution  as  x  •+  +  00 ,  and  thus  -  0  and 
k2  >  S2(l  -  ?2)/h3. 

Next,  in  Region  II,  h(x)  ■  h2  -  ^1  -  h2jx  such  that  6  =  h2  <  1 
and  y  *  -£l  —  h2j  <  0.  Then 

*  exp^-C2/2J  jA2M^a2,l,+C2j  +  B2U  |a2,l,+£2jj 


with  £ 2  ”  +2k|^fi2/£l-h2j  j  "  X]  ancl  a2  m  \  ~  \ 


-S32(l-f2)/2jl-E2jk, 


Finally,  in  Region  III,  h(x)  =  h2  +  ^  -  K2j (A/B)x  such  that 
6  -  h2  <  1  and  y  »  [^3 ”^2)  > 

Then  Cm  -  exP [+^3 ^ 2J “(A3M (a3 * 1  * ~ ^3)  +  B3U (a3 * 1 » "^3) } 

with  C3  =  -2kpi2b/pl3  -  h2j  +  xj  and  =  (l+f)/2-u>2  (l-f2)b/2kpi3  -  h2 


Continuity  of  free  surface  height  and  mass  flux  at  x  =  -1,  x  8  B/A  =  b 
requires 


I  <-« 

= 

I  <-« 

= 

Un(-D 

II  <°> 

= 

W°> 

II (0) 

St 

GIII(°) 

III(E> 

= 

h* (5) 

m(b) 

= 

uIV  (S) 

Using  Ul  =  k/l-h2,  p2  =  h2Pl,  u3  =  kh2b/h3-h2,  u4  =  kh3b/h3-h2 
and  ox  =  |«2>1,2mjJ  ;  «2  =  Ja2,l,2ji2J;  ^  =  Ja',1,2^}; 

62  =  ^a3,l,2u4j  and  M'=  dM/d£,  U'  =  dU/d£,  the  non-zero  values  of 

Cu  -  -expf-Aj 

C12  ‘  *“P(-Ul)K(ai) 

C13  "  +e,p(',‘l)U(“l) 


are 


-Xlexp(-Xl) 


C22=  k  expf-yj  [MfaJ 
C23=  kexpj-yjfujaj 

C32  =  exp(_1J2)M(a2) 

C33  =  exp(-P2)u(a2) 

C34  =  -exp[-w3]M(ei) 

C35  =  ~exp(_ll3)U(ei) 

C42  k  exp(-vi2)  [M(a2] 

C43  =  S  exp(-p2)[U(a2) 

C44  -+t  exp  (-1.3)^) 

C45  =  +k  exp(-p3)  [ufej 
C54  =  expf-y4)M(62] 

C55  =  exp|-u4ju|82j 

C56  =  "exp('X4S) 

C64  texPN[M(P2) 
C65  1  exp(-p4][U(62) 

C66  ~A4eXp (~A4E) 


m'N] 

2U'N] 


2”'N] 

2m'(bi)1 

20-  (Bjj 


2M'(B2)] 
20’ («2)] 


Again,  it  is  instructive  to  examine  the  case  of  no  rotation, 

f  *  0,  and  a  symmetric  topography,  A  =  B  and  =  H^.  The  cutoff 

points  along  the  boundary  curve  w  *  k  may  again  be  found  using 

X,  *  A,  »  0.  For  this  case,  a_  =  a'  =  a. 

14  2  3 


J.  V  H 

CIX  =  expj^-k2/^l-h2jj|A2M^a,l,2kz/^l-h2jj  +  B2U^a,l,2kz/ ^l-h2j  j  | 
?in  *  exp^-kz/^l-h2jj|A3M^a,l,2kz/^l-h2jj  +  B3U|a,l,2kz/ ^l-h2j 


where  z  =  -  ^l-h^jx  in  Region  II  and  +  ^l-h^jx  in  Region  III. 

The  continuity  conditions  require 


if  (s  =  "1}  " 0 


iP  (X  =  +d  - 0 


CrI  (x  =  0)  =  CZII  (x  =  0) 


d^n  d^III 

(,  =  °)  =  _I  (5E  ,  0) 


-M(ct) 


-U(a) 


M(a)-2Mf  (a)  U(a)-2U*(a)  M(a)-2Mf(a)  U(a)-2U' (a)  B„ 


M(B)-2M'(B)  U(B)-2Ui (B) 


M(B)-2M'(B)  U(B)-2Ut(B)  |B. 


where  a  =  «2  and  B  =  B2» 

The  determinant  of  these  coefficient s,  which  must  vanish  for  a 
nontrivial  solution,  may  be  reduced  to  a  product 

M(a) -2M* (a)  U(cO-2U'(a)  M(a)  U(a)  = 

M(B)-2M,(B)  U(B)-2Uf(B)  M(B)-2M1 (B)  U(B)-2U'(B) 


which  corresponds  to  symmetric  modes  (d£/dx  =0  at  x  =  0,1)  and 
antisymmetric  modes  (£  =  0  at  x  =  0,  d£/dx  =0  at  x  =  1) 
respectively. 

Numerical  results  are  shown  in  Figures  2-5  for  the 
symmetric  case  =  3  Km,  =  1  Km,  =  3  Km  and  A  =  B  =  60  Km  and 
a  slightly  asymmetric  case,  ^  *  3  Km,  H2  -  1  Km,  H3  =  4.85  Km  and 
A  *  60  Km,  B  =  78  Km,  both  without  rotation,  f  *  0.  Figure  2  shows 
the  symmetric  case  dispersion  equation  followed  by  Figure  3 
which  show  the  first  two  mode  shapes.  Note  that  both  symmetric  and 
antisymmetric  modes  result.  Figure  4  shows  tne  asymmetric  case 


dispersion  equation  followed  by  Figure  5  which  shows  the  first  two  mode 
shapes;  note  the  slight  shift  in  mode  shape  caused  by  the  asymmetry. 


COMPARISON  OF  CUTOFF  FREQUENCIES 

The  cutoff  frequencies  for  a  discontinuous,  segmented  constant 
depth  ridge  defined  by  depth  H  for  x  <  -B  ,  H  for  -B  <  x  <  B  and 

1  4  L  C 

for  x  >  +BC  are  given  by 

=  (mr/2)  Jh2/^l-h2jJ  1/2 

where  the  depth  and  half  ridge  width  are  again  used  as  reference 
lengths,  e.g.  Buchwald  (1968),  Shaw  and  Neu  (1979).  Corresponding  cutoff 
frequencies  may  be  found  for  the  linear  ridge  for  various  values  of  h2; 
these  are  shown  in  Table  Is 


TABLE  1  -  LINEAR  RIDGE  CUTOFF  FREQUENCIES,  wn|h2j 


hjsr 

1 

2 

3 

4 

5 

6 

0.1 

1.395 

3.220 

4.570 

6.189 

7.595 

9.166 

0.2 

1.731 

3.796 

5.502 

7.360 

9.094 

10.920 

0.3 

2.171 

4.394 

6.424 

8.541 

10.583 

12.678 

0.4 

2.580 

5.064 

7.438 

9.849 

12.223 

14.620 

0.5 

3.053 

5.862 

8.635 

11.400 

14.160 

0.6 

3.645 

6.881 

10.150 

13.368 

0.7 

4.456 

8.297 

12.249 

0.8 

5.800 

10.567 

The  ratio  of  these  cutoff  frequencies  to  those  for  the  constant 

depth  case  varies  with  h2  and  slightly  with  n  and  is  shown  in  Figure  6. 

Here,  u>n  (constant  depth) (linear  depth)  is  shown  as  a  function  of  h2 

with  a  range  indicated  at  each  h2  representing  the  variation  with  n. 

For  a  wide  range  of  values  of  h2,  the  ratio  is  approximately  0.5. 

This  implies  that  a  linear  ridge  of  total  width  2B  and  a  constant  depth 

L 

segment  ridge  of  one  half  of  this  width,  Bc  ■  0.5  B^,  would  lead  to 
approximately  the  same  cutoff  frequencies.  This  corresponds  to  maintaining 
the  same  total  cross  sectional  area  of  the  ridge  and  is  physically 
appealing  as  a  useful  approximation  allowing  a  complicated  problem  (linear 
ridge)  to  be  replaced  by  a  much  simpler  problem  (constant  depth  ridge). 


Of  course,  this  example  does  not  include  rotation  and  is  only  for  the 
symmetric  ridge.  Furthermore,  there  are  other  features  of  the  dispersion 
equation  which  should  also  be  matched,  e.g.  the  minimum  in  the  group 
velocity  noted  by  Buchwald  (1968).  However  the  initial  point  (cutoff 
frequency)  and  the  asymptote  (ai2  -  h2k2)  are  wel1  matched  by  this 
substitution.  This  is  particularly  gratifying  in  view  of  the  drastic 
difference  in  bottom  slopes  which  arise  in  these  two  different  cases. 


ACKNOWLEDGEMENT:  This  research  was  supported  by  the  Office  of  Naval 

Research  -  Physical  Oceanography  -  under  contract  N00014-79-C-0067 . 


kj 


REFERENCES 


BUCHWALD,  V.T.  1968a.  Long  waves  on  oceanic  ridges.  Proc.  R.  Soc.  Lond., 

A,  308;  343-354. 

ECKART,  C.  1951.  Surface  waves  in  water  of  variable  depth.  Marine  Physical 
Laboratory  of  the  Scripps  Institute  of  Oceanography,  Wave  Rep.  No.  100, 
S.I.O.  Ref.  51-12,  99  pp.  (unpubl.  manuscr.). 

GUZA,  R.T.  and  R.E.  DAVIS.  1974.  Excitation  of  edge  waves  by  waves  incident 
on  a  beach.  J.  Geophys.  Res.,  79:  1285-1291. 

HEATH,  R.  1979.  Tsunami  Amplification  in  New  Zealand  due  to  Trapped  Waves. 
IUGG  Tsunami  Symposium,  Canberra,  Dec.  1979  (to  be  published  in 
proceedings) . 

HIDAKA,  K.  1976.  Seiches  due  to  a  submarine  bank.  Topics  in  Ocean 

Engineering  III,  ed.  C.  Bretschneider ,  Gulf  Publishing  Co.,  Houston, 

Tx.,  pp.  137-140. 

LeBLOND,  P.H.  and  L.A.  MYSAK.  1978.  Waves  in  the  Ocean,  Elsevier 
Oceanography  Series,  Elsevier  Scientific  Publishing  Co.,  N.Y. 

LONGUET-HIGGINS,  M.S.  1969c.  On  the  trapping  of  long-period  waves  round 
islands.  J.  Fluid  Mech.,  37:  773-784. 

MEYER,  R.E.  1971.  Resonance  of  unbounded  water  bodies.  Mathematical 

Problems  in  the  Geophysical  Sciences,  ed.  W.H.  Reid,  Amer.  Math.  Soc., 
189-228. 


MONK,  W.H.,  F.E.  SNODGRASS  and  G.  CARRIER.  1956.  Edge  waves  on  the 
continental  shelf.  Science,  123;  127-132. 


i 


NEU,  W.  and  R.P.  SHAW.  1978.  Long  waves  on  linearly  varying  bottom 
topography.  Int.  Sym.  Long  Waves  in  Oceans,  188-192;  1979. 

ISLWO  Manuscript  Report  53,  Marine  Sciences  Directorate,  Ottawa, 
Canada,  188-192. 

REID,  R.O.  1958.  Effect  of  Coriolis  force  on  edge  waves,  I. 

Investigation  of  normal  modes.  J.  Mar.  Res.,  16:  109-144. 

ROBINSON,  A.R,  1964.  Continental  shelf  waves  and  the  response  of  sea 
level  to  weather  systems.  J.  Geophys.  Res.,  69:  367-368. 

SHAW,  R.P.  1974.  Long  waves  on  linear  topographies.  JTRE  Internal  Report 
No.  119,  Haw.  Inst,  of  Geophysics,  Honolulu,  Hi. 

SHAW,  R.P.  1977.  Long  waves  obliquely  incident  on  a  continental  slope  and 
shelf  with  a  partially  reflecting  coastline.  IUGG  Tsunami  Symposium, 
Ensenada,  Mexico  (also  see  1979  Marine  Geodesy  2(1):  1-14). 

SHAW,  R.P.  and  W.  NEU.  1979.  Long  wave  trapping  by  parallel  bottom 

contours.  Proceedings  of  IUGG  Tsunami  Symposium,  Canberra,  Australia 
(in  press) . 

SHEN,  M.C.,  R.E.  MEYER  AND  J.B.  KELLER.  1968.  Spectra  of  Water  Waves  in 
Channels  and  around  Islands.  Phy.  of  Fluids,  Vol.  11(11):  2289-2304. 

STOKES,  G.G.  1846.  Report  on  recent  researches  in  hydrodynamics. 

Rep.  16th  Meet.  Brit.  Assoc.  Adv.  Sci.,  Southampton,  1846. 

Murray,  London,  pp.  1-20.  (also  see  Math.  Phys.  Pap.,  1,  p,  167). 

URSELL,  F.  1952.  Edge  waves  on  a  sloping  beach.  Proc.  R.  Soc.  Lond.,  A, 
214:  79-97. 


TABLE  OF  FIGURES 


I 

I 

I 

I 


1.  Linear  Segment  Ridge  Approximation. 


2.  Symmetric  Peak — Dispersion  Equation  in  terras  of  the  physical  frequency 

o  in  rad/sec  and  k.  A  =  B  =  60,  =  3,  =  1  Km. 

3.  Symmetric  Peak — First  and  Second  Modes  for  various  values  of  k. 


4.  Asymmetric  Peak — Dispersion  Equation  in  terms  of  the  physical  frequency 
a  in  rad/sec  and  k.  A  =  60,  B  =  78,  =3,  H2  = 

H3  =  4.85  Km. 


5.  Asymmetric  Peak — First  and  Second  Modes  for  various  values  of  k. 


6.  Ratio  of  constant  depth  ridge  cutoff  frequencies  to  linear  ridge  cutoff 
frequencies  as  a  function  of  ridge  peak  to  base  ratio, 


Linear  Segment  Ridge  Approximation 


60-4-60 


J - 1 _ I _ I _ 1  I  l 

4  6  8  10  12  14  16 


k 


Figure  2 


Non- Symmetric  Peak 
Resonances 


SECURITY  CLASSIFICATION  OF  THIS  PA*'»F  (ith*n  Pntn  I’ntoriul) 


REPORT  CQCUMENT ATIOM  PAGE 


1  REPORT  NUMBER 


I  4  TITLE  (and  Subtitle) 


Hpr;  RliAD  INSTRUCTIONS 

_ DKFOKK  COMPUTHNO  KORM 

2.  GOVT  ACCESSION  NO.  *3  RECIPIENT'S  CATALOG  NUMBER 


1?0 


Long  Wave  Trapping  by  a  Linear  Ridge 


7.  author^; 


R.P.  Shaw  and  W.  Neu 


I  9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

Technical 

6  PERFORMING  ORG.  REPORT  NUMBER 

119  + 

8.  CONTRACT  OR  GRANT  NUMBERS 

N0001479  C0067  ^ 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 


Dept,  of  Eng.  Science 
SUNY  at  Buffalo 

Buffalo,  NY  14214  _ _ 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  REPORT  DATE 

Office  of  Naval  Research  -  Physical  Oceanography  August  1980 
Code  481  -  Bay  St.  Louis,  Mississippi  \ 3.  n um 0 e r  o f  p ages 

23 

Ti-'  MONITORING  AGENCY  NAME  A  ADDRESSflf  cUtforont  from  Controlling  Office)  15.  SECURITY  CLASS,  (of  thfa  report) 

Unc lassif ied 

15a.  dITcl  ASSIFICATION/DOWNGRAO  <n"g~ 

schedule 


13.  NUMBER  OF  PAGES 


I  16.  Ol  ST  HI  UU  Tl  ON  Si  ATEMENT  (of  this  Report) 


Approved  for  Public  Release 


Distribution  is  Unlimited 


17.  0ISTRI8HTI0N  STATEMENT  (of  the  ebatroct  entered  In  Block  20,  II  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WOROS  (Continue  on  rovoreo  aide  It  neceeeery  end  Identify  by  block  number) 


Edge  Waves 


20.  ABSTRACT  (Continue  on  reverae  aide  II  neceeeery  end  Identify  by  block  number) 

Long  waves  are  affected  by  bottom  topography  and  under  certain  conditions  may 
be  trapped  along  topographical  contours  which  then  act  as  wave  guides 
transmitting  wave  energy  for  great  distances  with  little  loss.  This  study 
examines  waves  trapped  along  a  submerged  ridge  described  by  straight  parallel 
bottom  contours  which  in  cross  section  are  composed  of  constant  slope  segments 
bounded  on  either  side  by  constant  depth  segments.  Solutions  are  found  for 
time  harmonic  waves  periodic  in  the  along  ridge  direction  and  of  exponential 


DO  ,'Srn  1473  EDITION  OF  1  NOV  65  1%  OOSOLETE 
$  N  0102-  LF*  0  I  4*  660  I 


SECURITY  CLASSIFICATION  OF  THIS  fWh*#t  Date  Knifed) 


decay  behavior  normal  to  the  ridge  over  the  constant  depth  segments. 

Over  the  linearly  varying  topography  describing  the  ridge,  the  solution  is 
in  terms  of  two  Kummer  (or  Whittaker)  functions.  For  a  given  geometry, 
a  dispersion  equation  is  obtained  relating  the  wave  frequency  to  the  along 
ridge  wave  number  for  trapped  waves.  A  constant  Coriolis  parameter  is 
included,  but  primary  interest  is  on  class  I  (high  frequency)  waves. 

A  comparison  of  cutoff  frequencies  predicted  for  this  piecewise  continuous 
ridge  and  those  for  a  segmented  constant  depth  ridge  is  made,  and  the 
appropriate  scaling  factors  between  the  two  results  are  discussed. 


